Persistent Current States in Bilayer Graphene 
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Broken symmetry states that carry equilibrium currents were proposed some time ago by re- 
searchers seeking a microscopic theory of superconductivity, but have never been observed. In 
this Letter we point out that the unusual electronic structure of unbalanced bilayer graphene is 
particularly favorable for the realization of this interesting type of ordered state. 



The search for a microscopic theory of superconduc- 
tivity, ultimately brought to a conclusion by the success- 
ful work of Bardeen, Cooper, and Schrieffer pQ, led a 
number of early condensed matter theory researchers [2]- 
[4] to speculate on the possibility that interactions could 
under some circumstances lead to electronic states that 
carry currents in equilibrium. Such a state would have 
to break rotational and time-reversal symmetry. From 
a more modern point of view, its continuous formation 
would be signaled by an angular momentum £ — 1 charge- 
channel Pomeranchuk [5 (Fermi-surface distortion) in- 
stability. 

An £ = 1 Fermi-surface distortion transfers occupa- 
tion between quasiparticle states with opposite current 
components along the continuously variable direction 
in which overall current flows, as illustrated schemati- 
cally in Fig. [T] For a translationally-invariant electronic 
system, however, it is known from both microscopic- 
quantum-mechanical |6j [7] and Fermi-liquid-theory [8] 
points of view that such a distortion corresponds sim- 
ply to a Galilean boost which raises the center-of-mass 
kinetic energy and does not change the interaction en- 
ergy. An £ = 1 charge-channel Pomeranchuk instability 
is therefore an impossibility in an electron fluid. In a 
crystal, where the Pomeranchuk instability notion refers 
to Fermi-surface distortions which reduce lattice symme- 
tries, £ = 1 distortions have also been viewed [9 as ex- 
tremely unlikely. Instabilities with £ = 2, which lead to 
electron nematic states, do on the other hand appear to 
occur [TQl [11] in a variety of different systems and spin- 
channel £ = instabilities, which lead to ferromagnetism, 
are common. In this Letter we point out that a Bernal- 
stacked graphene bilayer (BLG) in which a gap has been 
opened by a transverse electric field is ideally suited to 
host an £ = 1 Pomeranchuk instability because its elec- 
tron and hole Fermi seas are disks that are spread widely 
over momentum space, as illustrated in Fig. [TJb) [T2] . 
Fermi-surface distortions with £ = 1 can lower energy by 
compactifying the Fermi sea to realize the momentum- 
space condensation envisaged by London [4 more than 
seventy years ago. 

Recent progress [13] has made it possible to prepare 
and study the electronic properties of two-dimensional 
(2D) electron systems based on single- and few-layer 
graphene. This advance has provided researchers with 




FIG. 1: (color online) £ = 1 Pomeranchuk instability of a con- 
ventional Fermi liquid and of density- unbalanced BLG. The 
lower panels illustrate the quasiparticle band dispersions that 
correspond to the Fermi surfaces illustrated in the top pan- 
els, a) The circular orange Fermi circle with Fermi radius kpo 
denotes the occupied states of a two-dimensional Fermi sea. 
A Pomeranchuk instability occurs when energy is reduced by 
a momentum-space angle-dependent change in Fermi radius 
that is proportional to cos (-£#). For £ — 1 the distorted Fermi 
surface, indicated here by the dashed black circle, is simply 
shifted in momentum space, b) The disk-shaped orange re- 
gion with inner and outer Fermi radii /cFo,in and /cFo,out de- 
notes occupied conduction band states in unbalanced BLG. 
Shifts in the inner and outer Fermi radii that are proportional 
to cos(#) and have opposite signs for the inner and outer Fermi 
lines, thicken and thin the disk in opposite directions as indi- 
cated by the black dashed lines. The distorted state will tend 
to have lower interaction energy because the occupied states 
are more concentrated in momentum space. 

a new family of materials whose electronic structure is at 
the same time remarkably simple and remarkably vari- 
able. Single-layer graphene is described by a massless- 
Dirac-fermion model with conduction and valence bands 
that touch at two different points in momentum space 
and disperse linearly over a wide energy region. In a 
neutral Bernal-stacked bilayer with a very strong electric 
field directed perpendicular to the layers [14], the Fermi 
level lies within the conduction band of one layer and 
the valence band of the other. Inter-layer hybridization 
with strength 71 opens up an avoided crossing gap cen- 
tered at a finite 2D wavevector magnitude between states 
localized in opposite layers, yielding an unusual semicon- 
ductor with an electrically tunable gap. Our interest here 
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is in the electronic properties of degenerate electrons in 
the conduction band (or holes in the valence band) of 
this density-unbalanced configuration of BLG. When an- 
gular variation of the avoided crossing gap due to trigo- 
nal warping effects is neglected [15 j, the conduction band 
minimum occurs along a circle in momentum space and 
the non-interacting electron Fermi surface is an annulus 
as indicated in Fig. [TJb). Because the band density-of- 
states diverges as energy approaches the conduction band 
minimum it is clear, as noted previously p~6j [17], that 
interactions may play a central role in determining elec- 
tronic properties and that instabilities that break symme- 
tries are likely. The possibilities include common broken 
symmetries which lead to ferromagnetic or density-wave 
states but also, as we argue below, the spontaneous cur- 
rent states imagined by London, Heisenberg and others. 

We illustrate our main point by considering a toy 
model which ignores trigonal warping and spin and val- 
ley degrees-of-freedom, and by using mean-field theory 
to estimate the toy model Fermi-liquid parameters. The 
Hartree-Fock energy functional [18] is 
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Eu F [{n k }} = ^[eb(k) - v]n k - — ^ n h V h - W n h > , 

(i) 

where £h(k) are bare-band energies, which are isotropic 
and thus depend only on k = |fe|, and A is the 2D elec- 
tron system area. In our toy model we assume that the 
interaction V k - k > depends only on \k — k f \. Expanding 
this energy functional in powers of the deviation Sn k from 
the occupation numbers corresponding to the undis- 
torted Fermi surface yields an energy expression of the 
Fermi-liquid-theory form: 

EHF[{nk}} = Eq+^e^ ¥ (k)8n k -^- ^8n k V k - k >8n k > , 
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where the Hartree-Fock band energy £^ F (k) ^ s defined by 
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Because the total energy is sensitive to cancella- 
tions between quasiparticle-velocity renormalizations 
and quasiparticle-interaction corrections, we must treat 
the two contributions on an equal footing. Below we con- 
sider the quasiparticle-velocity renormalization first and 
use this analysis to define our notation. 

We begin by calculating the Hartree-Fock velocities, 
i.e. the momentum-space radial derivatives of the quasi- 
particle energy at the outer (n — out) and inner (n = in) 
Fermi circles. The velocities are most conveniently eval- 
uated by replacing the derivative with respect to k in 
V k - k > by a derivative with respect to k' . When the in- 
teraction correction to the velocity is integrated by parts, 
this derivative then acts on n^) picking out states at the 



Fermi energy. We find that 

v™(k) = sgn(n)K(fc)| + ^ j J cos(^) 

x {[fc / V r fe_fc'] fe ' = fe FOout - [k'V k -k'] k > = k FO in } , 

(4) 

where (j) is the difference between the angular coordinates 
of k and k' and we have noted that V k - k > depends only 
on (j) once k and k' are fixed [19]. This leads to 
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where v* ut = ^(fc F o,out), "out = e 2 /(Hv out ), uffl = 

Vm(kYQ,n, &F0,n') (&F0,n k F0 , n ') 1/2 /( 27re2 ) is a dimension- 
less interaction parameter that is symmetric in the n, n' 
inner/outer indices, and the Kn's are Fourier components 
of the interaction's (^-dependence. The corresponding ex- 
pression for v* n can be obtained by interchanging the in 
and out labels. 

To look for Pomeranchuk instabilities we paramaterize 
the inner (n = in) and outer (n = out) Fermi surfaces in 
terms of dimensionless distortion functions: 

&F,n = &F0,n[l + CL n (0)] = fc F 0,n + £&F,n(#) , (6) 

and expand these in terms of their angular momentum 
components a n (6) = Yl~i^-oo a n£e t£e . The distorted state 
has a 5n k which is non-zero only in the vicinity of the 
inner and outer Fermi lines. It is easy to verify that the 
first-order correction to the energy vanishes. To obtain 
the second-order correction we linearize the Hartree-Fock 
energies around &F0,in and &F0,out : 



~hv* n (k - fc F0 ,in) 
ft V out( k ~ ^F0,out) 



and add the quasiparticle interaction contribution to ob- 
tain the energy change for small distortions: 



h 



+oo 

-T y 

47T ^ ^ 

n.n' £— — oo 



-sg: 



i n (nn')C'} 4o?n» *Efl?n' o„fa,V< 



(7) 



This is our principal result. Note that for a Galilean- 
invariant system which has only an outer Fermi radius 
the two interaction contributions to the £ = 1 distortion 
energy cancel, recovering the general result [6HH] men- 
tioned previously. An £ = 1 Pomeranchuk instability 
occurs when the determinant of the following matrix 
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FIG. 2: (color online) The dimensionless "outer-inner" pseu- 
dopotential U^ t in vs. e = b/ (2k) for £ = (s-wave) and ^ = 1 
(p-wave). This plot is for 5 = 2/(a e fik) = 1 where a e ff and /c 
are defined in the main text. The parameter b measures the 
thickness of the ring and b — >• in the limit in which the Fermi 
energy sf equals the band-edge energy e m i n . The left panel 
is for a very small screening parameter, A = 0.01, while the 
right panel is for fully screened Thomas- Fermi interactions, 
A = 1.0. 



is zero. Using Eq. ([5| (and the corresponding equation 
for v* n ) we finally find the following criterion: 
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Notice that this instability criterion depends only on the 
"outer-inner" interaction. 

We estimate the disk Fermi-liquid interaction parame- 
ters using 
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where <2tf is the Thomas-Fermi screening wave vector, 
and A G [0, 1] is a dimensionless control parameter that 
allows us to interpolate between bare-Coulomb (A = 0) 
and Thomas- Fermi (A = 1) limits. The Thomas- Fermi 
screening wave vector is proportional to the density-of- 
states at the Fermi energv, N(0) = 2m e fik/(7rH 2 b) where 
b = &F0,out - &F0,in, k = (&F0,out +&F0,in)/2, and m e ff pa- 
rameterize the band-energy dispersion at its minimum. 
It follows that g T F = 27re 2 7V(0) = 4k/(a e ^b) where 
a e ff = h 2 /(m e ffe 2 ) is an effective Bohr radius. The inner- 
outer interaction parameters U^ t in obtained using this 
approximation depend on two dimensionless quantities, 
e = b/ (2k) and 5 = 2/(a e fffc), and are plotted for £ = 0, 1 
in Fig. |2j The strong declines for small e are due to 
strong Thomas-Fermi screening and are certainly an ar- 
tifact of this approximation in combination with the di- 
vergent Fermi-level density-of- states for e — >• 0. We can 
make two conclusions based on these estimates: i) Inter- 
actions in the i = 1 channel are strong enough to produce 
an £ = 1 Pomeranchuk instability and ii) Interactions in 
the £ — channel are likely even stronger than those in 
the £ = 1 channel. If bilayer graphene had a single val- 
ley, £ = instabilities, which imply phase separation into 
low and high density regions, would be forbidden by the 
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FIG. 3: (color online) Quasiparticle bands, Hall conductivi- 
ties and spontaneous currents due to £ = and £ = 1 Pomer- 
anchuk instabilities in unbalanced BLG. These results were 
obtained with an external potential difference corresponding 
to a perpendicular electric field of 1 V/nm between layers, (a) 
Valley polarized quasiparticle bands for three different total 
carrier densities (both spins). The upper and lower panels il- 
lustrate conduction band energy dispersions and constant en- 
ergy contours respectively. The left and right columns are for 
momenta near the K' and K valleys. The £ — instability is 
signaled by unequal occupation of the two valleys. The £ — 1 
instability is signaled by the breaking of rotational symmetry 
within an occupied valley. The solid black line marks the po- 
sition of the Fermi level which is chosen as the zero of energy: 
£f = 0. (b) Spontaneous Hall effect vs. carrier density. The 
Hall effect [Eq. ( |11[ )] reflects spontaneous valley polarization. 
We find that the Hall conductivity increases monotonically 
with carrier density until minority (K') valley with opposite 
Hall conductivity starts to be occupied. The shaded regions 
indicate intervals of carrier density over which spontaneous 
currents occur due to £ — 1 Pomeranchuk instabilities. The 
lower density instability occurs within the majority valley and 
the higher density instability within the minority valley, (c) 
Total spontaneous persistent current evaluated from Eq. (12) 
for both spins. The conduction and valence bands of the po- 
larized valley have similar contributions to the total persistent 
current. For this calculation we used a /c-point sampling den- 
sity near the Dirac cones equivalent to a density of 4608 x 4608 
points in the whole Brillouin zone. The blue dotted line rep- 
resents a guide to the eye for the numerical data points. The 
slope of the continuous line indicates the effective quasipar- 
ticle velocity of the current carrying states, which is about 
1.5% of the Fermi velocity in graphene. 



long-range Coulomb interaction. However, because of the 
presence of valley and spin degrees of freedom it appears 
that the strongest instabilities of the normal state will 
have £ — 0, leading to spin [16 or valley polarized states. 
£ = 1 Pomeranchuk instabilities do however appear pos- 
sible within these polarized states. 

We have corroborated the conclusions drawn from our 
toy-model calculation by performing self-consistent tt- 
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orbital lattice Hartree-Fock calculations similar to those 
described in Ref. [20j For these calculations we added in- 
teractions to a band model with nearest-neighbor intra- 
layer hopping 70 = 3.12 eV and inter-layer hopping 
71 = 0.377 eV, ignoring other hopping parameters for 
the sake of simplicity [21]. We have suppressed spin- 
polarization instabilities in these calculation so that i = 
instabilities are manifested by valley rather than spin- 
polarization. Because external electric fields in unbal- 
anced BLG lead to large Berry curvatures fl nj k of oppo- 
site sign in the vicinity of K and K' valley points [22] 
the anomalous Hall conductivity can be used as an ob- 
servable for spontaneous valley polarization. The Hall 
conductivity (per spin) is calculated by integrating ft Uy k 
over occupied quasiparticle states near the Dirac points, 

a H = -r I — —2 ^ /n,fc ^n,fc , (H) 
n J ( 27 n 

where f n ^ is the Fermi-Dirac distribution and a sum is 
carried over the band index n. 

To evaluate the spontaneous current density associated 
with the t = 1 instability we use the current density oper- 
ator j = — eVkHo/(AK). Because the interaction Hamil- 
tonian is independent of position, the current operator 
is given by its non-interacting one-body form, which is 
diagonal in Bloch momentum. Its expectation value (per 
spin) in the mean-field state is therefore 

(j) = -Jj^^2fn,k (^n,fc|V fc H |^n,fc) , (12) 
k , n 

where \(p n ,k) are the Hartree-Fock single-particle states. 
We evaluate the current density averaging the numerical 
derivative of the tight-binding Hamiltonian in momen- 
tum space. Fig. [3] plots Hall conductivities and sponta- 
neous currents evaluated for this model which demon- 
strate that spontaneous valley polarization occurs for 
carrier densities smaller than ~ 25 x 10 11 cm -2 , and 
spontaneous currents for carrier densities smaller than 
~ 9 x 10 11 cm -2 and again near the onset of minority 
valley occupation. These mean-field calculations likely 
overestimate the density-ranges over which the broken 
symmetry states occur. 

Because they appear only at low carrier densities, 
comparable to or smaller than typical disorder-induced 
density-fluctuation scales for bilayer samples on silicon 
oxide substrates, and because large electric fields are fa- 
vorable for their occurrence, we anticipate that momen- 
tum space condensation is at present most likely to oc- 
cur in dual-gated bilayer graphene samples on h-BN sub- 
strates [23] . We expect that trigonal warping of the un- 
balanced bilayer conduction bands will favor momentum 
space condensation over competing [17] density-wave in- 
stabilities. Because the i = 1 Pomeranchuk instability 
is likely only within states in which spin or valley po- 
larization, or both, has already occurred, its appearance 



should be signaled most clearly by observables that re- 
flect broken inversion symmetry, rather than simply bro- 
ken time reversal symmetry. Since spontaneous currents 
cannot flow in closed samples, or even in samples that are 
in contact with an unbiased external circuit, equilibrium 
momentum space condensation must be accompanied by 
domain formation. We expect that this property will 
lead, not to dissipationless current flow of the type ob- 
served in superconducting samples, but instead to highly 
non-linear and hysteretic current- volt age characteristics. 

Heisenberg and London imagined that momentum- 
space condensation might explain superconductivity. If 
spontaneous current states do occur in unbalanced bi- 
layer graphene, it will be fascinating to sort out the sim- 
ilarities and differences between the properties of these 
broken symmetry states and those of superconductors. 
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